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A simple confining separable interaction Ansatz for the rainbow-ladder truncated 
QCD Dyson-Schwinger equations is used to implement chiral restoration and quark 
deconfinement in a study of qq meson states at finite temperature. The model is 
fixed at T = by reproducing selected it and p properties. Deconfinement and 
chiral restoration are found to both occur at Tc = 146 MeV. In the pion sector, 
we investigate Mt^{T) and f-wiT) along with the exact QCD mass relation and 
the GMOR relation. For the vector mode, we investigate the 3-space transverse 
and longitudinal masses Mj(T) and M^{T), along with the width for the decay 
p'^ — ► e+e~. The equation of state (EOS) for the model is investigated in the 
T — fi plane, a tri-critical point is identified and the relationship to a bag model 
is discussed. The deconfinement transition in rapidly rotating neutron stars is 
considered and a new signal from the pulsar timing in binary systems with mass 
accretion is suggested. The model is used to discuss the possibility of a supercon- 
ducting quark matter phase. 

1 Introduction 

The experimental search for the QCD deconfinement phase transition in ul- 
trarelativistic heavy-ion coUisions wiU enter a new stage when the relativistic 
heavy-ion coUider (RHIC) at Brookhaven wiU provide data complementary 
to those from the CERN SPS El. It is desirable to have a continuum field- 
theoretical modeling of quark deconfinement and chiral restoration at finite 
temperature and density (or chemical potential /i) that can be extended also 
to hadronic observables in a rapid and transparent way. Significant steps in this 
direction have recently been taken through a continuum approach to QCD^^^ 
based on the truncated Dyson-Schwinger eauations (DSEs) within the Mat- 
subara formalism EIu'lI and a recent review a is available. A most appealing 
feature of this approach to modeling nonperturbative QCD^^^t is that dynam- 
ical chiral symmetry breaking and confinement is embodied in the the model 
gluon 2-point function constrained by chiral observables at T = /.t = and no 
new parameters are needed for extension to T, /i > 0. Approximations intro- 
duced by a specific truncation scheme for the set of DSEs can be systematically 
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relaxed. However due to the breaking of 0(4) symmetry and the number of 
discrete Matsubara modes needed, the finite T, /i extension of reahstic DSE 
models entails solution of a complicated set of coupled integral equations. The 
generation of hadronic observables from such solutions, although a straight- 
forward adaption of the approachQ found to be successful at T = /i = 0, adds 
further to the difficulties. In the separable model we study here, detailed real- 
ism is sacrificed in the hope that the dominant and essential features may be 
captured in a simple and transparent format. To this end we simplifiy an exist- 
ing T = /i = confining separable interaction Ansatz u to produce a gaussian 
separable model for T, /i > 

2 Confining separable Dyson- Sciiwinger equation model 

In a Feynman-like gauge where we take = 5^i,D{p — q) to be the effective 
interaction between quark colored vector currents, the rainbow approximation 
to the DSE for the quark propagator S{p) — [i'^A{p) + B{p) + mo]~^ yields in 
Euclidean metric 

Bip) ^ llf^Dip^q) m±I!^ (1) 

3 7(271)4 ^'q^A^{q) + [B{q)+mo]' ' ^' 

[Aip)^l]p' = If-pLoip-q) ^l^Mi^^ (2) 

^ 3 7 (2^)4 ^'q^A^{q) + [B{q)+mof 

We study a separable interaction given by0 

D{p-q) = Do k{p^)h{q^) + D, h{p^){p ■ q)h{q^) , (3) 

where Do,Di are strength parameters and the form factors, for simplicity, are 
here taken to be fi{p'^) = exp(— p^/A|) with range parameters A^. It is easily 
verified that if Dq is non-zero, then B{p) — Am /o(p^), and if Di is non-zero, 
then A{p) = 1 + Aa /i(p^). The DSE then reduces to nonlinear equations for 
the constants Am and Aa. The form factors should be chosen to simulate the 
p^ dependence of A(p) and B(ji) from a more realistic interaction. We restrict 
our considerations here to the rank-1 case where Di — and A{p) = 1. The 
parameters Dq, Aq and mp are used to produce reasonable tt and lo properties 
as well as to ensure the produced Blp) has a reasonable strength with a range 
Ao - 0.6 . . . 0.8 GeV to be realistic Q. 

If there are no solutions to p'^A'^{p) + (B{p) + mp)^ = for real p^ then 
the quarks are confined. If in the chiral limit (mp — 0) there is a nontrivial 
solution for B{p), then chiral symmetry is dynamical broken. Both phenomena 
can be implemented in this separable model. In the chiral limit, the model is 
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confining if Do is strong enough to make Am/Ao > l/\/2c- Tlius for a typical 
range Aq, confinement will typically occur with M{p ~ 0) > 300 MeV. 

Mesons as qq bound states are described by the Bethe-Salpeter equation 
which in the ladder approximation for the present approach is 

- A(p2)r(p, P)^lJ -0^Dip - q)^,Siq+)Tiq, F)5(«_)7^ , (4) 

where q± — q ± P/2 and P is the meson momentum. The meson mass is 
identified from A(P^ = —JVP) =1. With the rank-ljSeparable interaction, only 
the 75 and the 75/* covariants contribute to the ttLI, and here we retain only 
the dominant term T.^{p,P) = i'j^ETrip, P). For the vector meson, the only 
surviving form is Vp^{p, P) = 'y'^{P)Ep{p, P), with 'jJj^iP) being the projection 
of 7^ transverse to P. The separable solutions have the form Ei {p, P) = 
/o(p^)Ci(P^), i — TT,p, where the Ci factor out from Eq. (^). 

In the limit where a zero momentum range for the interaction is simulated 
by /o(<Z^) fX(5''(g), then the expressions for the various BSE eigenvalues A(P^) 
reduce to those of the Munczek and Nemirovskycl model which implements 
extreme infrared dominance via D{p — q) cx q). The correspondence 

is not complete because the quark DSE solution in this model has A{p) ^ 1. 
The T, /i > generalization of this infrared dominant (ID) model have been 
studied recently Ero. 

3 Pion and rho-meson properties 

With parameters mo/Ao = 0.0096, DqAI = 128 and Aq = 0.687 GeV, the 
present Gaussian separable (GSM) model yields M^r — 0.14 GeV, Mp — ~ 
0.783 GeV, U = 0.104 GeV, a chiral quark condensate (qq)^/^ = -0.248 GeV, 
and a p — 7 coupling constant gp ~ 5.04. 

The generalization to T 7^ is systematically accomplished by transcrip- 
tion of the Euclidean quark 4-momentum via q qn = {^n, <l), where w„ — {2n + l)7rT 
are the discrete Matsubara frequencies. The obtained T-dependence of the 
mass gap Am(T) allows for a study of the deconfinement and chiral restora- 
tion features of this model. We find that both occur at Tc = 146 MeV where, 
in the chiral limit, both Ato(T) and (qq)^ vanish sharply as (1 — T/Tc)^ with 
the critical exponent having the mean field value (3 = 1/2. 

For the qq meson modes, the 0(4) symmetry is broken and the type of 
mass shell condition employed must be specified. If there is a bound state, the 
associated pole contribution to the relevant qq propagator or current correlator 
will have a denominator proportional to 

1 - X{nl, P^) cx r!^ + P2 + M\T) . (5) 
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We investigate the meson mode eigenvalues A using only the lowest meson 
Matsubara mode (fim — 0) and the continuation — > The masses so 

identified are spatial screening masses corresponding to a behavior exp(— Mx) 
in the conjugate 3-space coordinate x and should correspond to the lowest 
bound state if one exists. 




Figure 1: T-dependence of qq meson properties for the rank-1 separable model up to 
Tc = 146 MeV. Left panel; Spatial Mt^{T) and quantities associated with mass relations as 
dictated by chiral symmetry; Right panel; Spatial masses Mj{T) and M^{T) along with 
the — > e+e~ decay width and the associated vector coupling constant gp{T). 



The obtained tt and p masses displayed in Fig. |l| and are seen to be only 
weakly T-dependent until near Tc ~ 146 MeV. This result for M-j^lT) repro- 
duces the similar behavior obtained from the ladder-xainbow truncation of the 
DSE-BSE complex with a more realistic interactiona. The qualitative behav- 
ior obtained for the 3-spaca transverse and longitudinal masses Mj(T), M^(T) 
agrees with that reportedu for the limiting case of the zero momentum range 
or ID model. 

To explore the extent to which the model respects the detailed constraints 
frorn_chiral symmetry, we investigate the exact QCD pseudoscalar mass rela- 
tion Ei3 which, after extension to T > 0, is 

Ml[T)U{T)^2morp{T) . (6) 

Here rp is the residue at the pion pole in the pseudoscalar vertex, and in the 
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present model, is given by 



irp{T) = tr, J ^ 755(<z„ + ^)T^{q,,; P)S{q„ ~ |) . (7) 

The relation in Eq. is a consequence of the pion pole structure of the isovec- 
tor axial Ward identity which linkiL the quark propagator, the pseudoscalar 
vertex and the axial vector vertex E2l. In the chiral limit, rp (qq)^ / and 
Eq.(^, for small mass, produces the Gell-Mann-Oakes-Renner (GMOR) rela- 
tion. 

The exact mass relation, Eq. (^) , can only be approximately satisfied when 
the various quantities are obtained approximately such as in the present sep- 
arable model. The error can be used to assess the reliability of the present 
approach to modeling the behavior of the pseudoscalar bound state as the tem- 
perature is raised towards Tc. Our findings are displayed in Fig. [|. There the 
solid line represents rp(T) calculated from the quark loop integral in Eq. (7); 
the dotted line represents rp extracted from the other quantities in Eq. (3). 
It is surprising that the separable model obeys this exact QCD mass relation 
to better than 1% for the complete temperature range up to the restoration of 
chiral symmetry. Also evident from Fig. ^is that M.^{T) and the chiral conden- 
sate are largely temperature independent until within about 0.8 Tc whereafter 

rises roughly as fast as the condensate falls. 

We have also investigated the (approximate) GMOR relation for the present 
model. The quantity (qq)^ /N^ is displayed in Fig. |l| as the long-dashed line, 
and if the GMOR relation were exactly obeyed, this would coincide with 

/7r/2m which is the dotted line. The quantity N'^ enters here via its role 
as the normalization constant of the chiral limit tt BS amplitude = 
ij5Bo{p'^)/N^. If all covariants for the pion were to be retaixifid and the axial 
vector Ward identity were obeyed, one would have N^^ = fir^- The results in 
Fig. indicate that the GMOR relation contains an error of about 5% when 
compared either to the exact mass relation or to the quantities produced by the 
separable model and that this is temperature-independent until about 0.9 Tc- 
It should be noted that f^,N'^ and are equivalent order parameters near 
Tc and have weak T-dependence below Tc- A consequence is that /^r, rp 
and {qq)'~' /N^ are almost T- independent and so are the estimated errors for 
the two mass relations linking these quantities. Since we obtain Af^ and /^r 
from the model BSE solutions at finite current quark mass, /^r does not exactly 
decrease to zero and Af^r does not exactly diverge at Tc- 

Vector mesons play an important role as precursors to di-lepton events 
in relativistic heavy-ion collisions and it is important to explore the intrinsic 
T-dependence of electromagnetic and leptonic vector coupling constants that 
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can arise from the quark-giuon dynamics that underHes the finite extent of 
the vector qq modes. The present model provides a simple framework for such 
investigations. The electromagnetic decay constant gpJ^) that describes the 
coupling of the transverse p° to the photon is given byEJ 



after accounting for the normalization i of the present BS amplitudes. The 
electromagnetic decay width of the transverse p mode is calculated from 



rpo^e+e-(r) = , . (9) 



At T = the experimental value is TpO_e+^- (0) — 6.77 keV corresponding to 
the value 5p(0) = 5.03. Our results for gp{T) and FpO^e+e-CT") are displayed 
in Fig. 1^. This electromagnetic width of the 3-space transverse p increases 
with T and reaches a maximum of 1.72 times the T = width at about 0.9 Tc- 
An increasing electromagnetic width for the p has been found empirically to 
be one oL-the possible medium effects that influence the heavy-ion dilepton 
spectrum 113. 



4 Equation of state (EOS) for quark matter 

The thermodynamical properties of the confining quark model and in particular 
the EOS and the phase diagram can be obtained from the grand canonical 
thermodynamical potential n{T,V,p) — T\nZ{T,V,p) = —p{T,p)V, where 
the contributions to the pressure (for a homogeneous system) 

p{T, p) = _Pcond(7', + Pkin{T, p) + po (10) 

are obtained from a mean-field approximation to the Euclidean path integral 
representation of the grand canonical partition function Z{T,V, p). In the 
rank-one separable gluon propagator model, the condensate contribution is 
Pcond{T,p) — 3 Am{T, p)^ /{16 Dq) and the kinetic part of the quark pressure 
is given by 

P„„,T,,0 = 2N,N, I Jfi^rpn H-P„..(T,., . (U, 

In Eq. kl = [{2n + \)tiT + ipf + k^, M{kl) = mo + Am{T,p)fo(kl) 

and the divergent 3-momentum integration has been regulaxized by subtract- 
ing the free quark pressure and adding it in the well-knownll3 regularized form 
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PficciT , fi) . The pressure contribution pq is found such that the total pressure 
( p!o| ) at the phase boundary in the T, /i-plane vanishes, see Fig. |. While in- 
vestigating this EOS for the separable confining quark model defined above we 
have observed that, as a function of the coupling parameter Dq, an instability 
d(pcond +Pkin)/dT < occurs whcn the criterion for confinement (absence 
of quasiparticle mass poles) is fulfilled. The physical quark pressure in the 
confinement domain of the phase diagram vanishes, see Fig. |^. The results 
for the EOS and the phase diagram can be compared to those for the zero 
momentum range model tl with the important modification that in the present 
finite range model the tricritical point is obtained at finite chemical potential 
whereas with zero range it was found on the fi = axis of the phase diagram, 
see Fig. ^ The location of the tricitical point could be experimentally veri- 
fied in CERN-SPS experiments provided that changes in the pion momentum 
correlation function could be detected as a function of the beam energyEJ. A 
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Figure 2: The quark matter pressure as a function of temperature for the separable confining 
model. The kinetic part (long dashed line) is overcompensated by the condensate part (dot- 
dashed line) resulting in an instability d(pt.ond +Pkin)/dT < (dotted line) for T < Tc = 
146 MeV which is characteristic for confining quark models. The physical pressure of quark 
matter (solid line) vanishes in this region. For comparison, the free quark matter pressure 
is shown by the dashed line. 

particularly interesting phenomenological application is the T — quark mat- 
ter EOS which is a necessary element for studies of quark deconfinement in 




7 



neutron stars. The present Gaussian separable model leads to a bag model 
EOS for quark matter with a bag constant B{T = 0) = 150 MeV/fm^ for the 
parameter set {Dq Aq = 128) employed in Sec. ^. A second parameter set 
{Do Aq — 97) that is also confining and provides an equally good description 
of the same tt and p properties produces B{T = 0) = 75 MeV/fm^; see also 
Fig. ^. More stringent constraints on the low-temperature EOS will require 
the inclusion of hadronic excitations including the nucleon. 
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Figure 3: Left panel: pressure (solid line) and energy density (long dashed line) vs. chemical 
potential at T = for the GSM with Aq = 0.756 GeV and DqA^ = 97.0 in the chiral 
limit. The phase transition from the confining quark matter phase to the deconfined one 
with restored chiral symmetry occurs at fi = 0.337 GeV and is first order. The results are 
coincident with a bag model EOS (dashed and dotted lines) for a bag constant B = 75 
MeV/fm'^. Right panel: quark matter phase diagram for the GSM with Aq = 0.687 GeV 
and DqAq = 128. Along the dashed line the pressure vanishes, the solid line separates the 
chiral symmetric phase from the broken one. A tricritical point is obtained at T ~ 127 MeV, 
/X ~ 120 MeV. 

5 Deconfinement in rotating neutron stars 

For the discussion of deconfinement in neutron stars, it is crucial to go beyond 
the mean field description and to include into the EOS also the hadronic bound 
states (neutrons, protons, mesons) in the confined phase. This task has not 
yet been solved aiiii therefore, we adopt for this phase a Walecka model as 
it is introduced inEj. In constructing the phase transition to the deconfined 
quark matter as described by the GSM with a parameter set leading to i? = 75 
MeV/fm"^ we have to obeypthe constraints of global baryon number conserva- 
tion and charge neutrality Il3. The composition of the neutron star matter is 
also constrained by the processes which establish (3— equilibrium in the quark 
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phase {d u + + D) and in the hadronic phase {n p + + D). For 
the given EOS we obtain a deconfinement transition of first order where the 
hadronic phase is separated from the quark matter one by a mixed phase in 
the density intervaLl.39 < n/riQ < 2.37, where uq = 0.16 fm^^ is the nuclear 
saturation densityllZl. 

Ah observable consequences should be discussed for fastly rotating com- 
pact objects and therefore we have studied these rotating configurations using 
this model-EOS with a deconfinement transition. The result is shown in Fig. 
^ and shows that within the present approach a deconfinement transition in 
compact stars is compatible with constraints on radius and mass recently de.- 
rived from the observation of QPO in low mass X-ray binaries (LMXBs) 
The basic quantity for the study of a deconfinement transition rotating com- 
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Figure 4: The mass M as a function of the equatorial radius Re (left panel) and the central 
density (right panel) for a neutron star with deconfinement transition. Rotating configura- 
tions (dashed lines) with maximum rotation frequency are connected with the static ones 
(solid lines) by lines of constant baryon number N/Nq = 1.3, 1.55, 1.8, 2.14, respectively. 
The occurence of an extended quark matter core in the compact star is compatible with re- 
cently derived constraints on maximum mass and radius from QPO observations in low-mass 
X-ray binaries, see text. 

pact stars is the moment of inertia which governs the rotation and thus the 
spin-down characteristics. Changes in the angular velocity Qlt) as a function 
of tixne can occur, e.g., due to magnetic dipole radiation t£l or mass accre- 
tionllJ. During the time evolution of an isolated pulsar, the deviation of the 
braking index from n{fl) = 3|-can signal not only the occurence of a quark 
matter core, but also its size t3. We have found that in LMXBs with mass 
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accretion at conserved angular momentum the occurence of a qupi matter 
core would reflect itself in a change from a spin-down to spin-up eratZI, see Fig. 
^. More detailed investigations of these scenarios have to be performed with a 



5 

a 3 

c 

1 



N./N„: 



i \ 2.14 
1 jj 1,85 



1400 .-7 


J lM„km'kHz]= _ 






eoo™«"°-^-- 




... 500 X \ 




400 














300 




00 




00 



n[kHz] 



1 1,2 1,4 1,6 1,8 2 22 2,4 
N./N„ 



Figure 5: Deconfinement signals from the pulse timing of rapidly rotating compact stars. 
Left panel: The braking index deviates from n = 3 when during the spin-down evolution 
of a pulsar a quark matter core is formed which entails a change in the moment of inertia. 
The larger the radius of the quark core, the more pronounced the signal. Right panel: The 
spin-down evolution of a compact star with mass accretion at constant angular momentum 
J flips to a spin-up behaviour at the onset of deconfinement. 

more realistic EOS, in particular for the hadronic phase. The possibility of a 
T = quark matter EOS which corresponds to a bag model EOS with small 
bag constants of the order of B = 70 MeV/fm'^ is an important result of the 
study of the confining quark model which bears interesting consequences for 
the study of further nontrivial phases in high-density QCD, as e.g. (color-) 
superconductivity. 

6 Superconducting quark matter 

The possible occurence of a superconducting quark matter phase has been 
recently reconsidered on the basiax£_D.onperturbative approaches to the effec- 
tive quark 4-point interaction Ei'LiiEjH and critical temperatures of the order 
of 50 MeV with diquark pairing gaps of w 100 MeV have been obtained. So, 
if quark matter occurs in the interior of compact stars as advocated in the 
previous section, then it had to be realised in such a superconducting phase. 
Deconfinement in compact stars can thus result in effects on the magnetic field 



10 



structjare as well as the cooling curves of pulsars. Contrary to previous esti- 
matesLJ, low temperature quark matter is a superconductor of second kind and 
thus the magnetic field can penetrate into the quark core in Abrikosov vortices 
and does not decay at timescales shorter than 10^ years. Thus the occurence 
of superconducting_quark matter phase in compact stars does not contradict 
observational datacj. The recently developed nonperturbative approaches to 
diquark condensates in high-density quark matter Edc3a can be further con- 
strained by studying the consequences for cooling_f;urves of pulsars c3 which 
have to be consistent with the observational datacEl. 
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A simple confining separable interaction Ansatz for the rainbow-ladder trun- 
cated QCD Dyson-Schwinger equations is found capable of modeling qq meson 
states at T > together with quark deconfinement and chiral restoration. De- 
confinement and chiral restoration are found to both occur at Tc — 146 MeV. 
The spatial screening masses for the meson modes are obtained. We find that, 
until near Tc, M^(T) and /^(T') are weakly T-dependent and that this model 
obeys the exact QCD pseudoscalar mass relation to better than 1%. The 
GMOR relation is found to be accurate to within 5% until very near Tc- For 
the vector mode, the 3-space transverse and longitudinal masses AlJ{T) and 
Mp{T) are weakly T-dependent while the width for the electromagnetic decay 
e+e~ is found to increase to 1.72 times the T = width. The equa- 
tion of state (EOS) for the model is investigated in the T — ji plane and it 
shows a tricritical point at T = 127 MeV, = 120 MeV. At T = the EOS 
can be given the form of a bag model where a broad range of bag constants 
B = 75 . . . 150 MeV/fm'^ is obtained consistent with possible parametrizations 
of vr and p observables. 

The consequences for deconfinement transition in rapidly rotating neutron 
stars are considered and a new signal from the pulsar timing in binary systems 
with mass accretion is suggested. The model EOS under consideration meets 
the new constraints for maximum mass and radius recently derived from QPO 
observations. Within the present model, quark matter below Tc ^ 50 MeV is 
a superconductor of second kind and it is suggested that the magnetic field 
in a neutron star forms an Abrikosov vortex lattice which penetrates into the 
quark matter core and thus in accordance with the observation does not decay 
on timescales of lO** years as previously suggested. 
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